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We consider the evolution of a dilute Bose-Einstein con- 
densate in an optical trap formed by a doughnut laser mode. 
By solving a one dimensional Gross-Pitaevskii equation and 
looking at the variance and the statistical entropy associated 
with the position of the system we can study the dynamical 
behavior of the system. It is shown that for small conden- 
sates nonlinear revivals of the macroscopic wave function are 
expected. For sufficiently large and dense condensates irre- 
versible dynamics takes place for which revivals of regular 
dynamics appear as predicted in JsjJ. These results are con- 
firmed by a two dimensional simulation in which the scales of 
energy associated with the two different directions mimic the 
experimental situation. 
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I. INTRODUCTION 

The experimental achievement of Bose-Einstein con- 
densation (BEC) in cold alkali atomic samples [[[|]3] has 
opened up new directions in atom optics. The develop- 
ment of optical (laser) devices to manipulate coherently 
atomic beams has been particularly spectacular. In the 
recent years it has led to the advances in atom cooling, 
providing thus the first steps towards the realization of 
coherent matter waves. Atomic beam splitters, or atomic 
mirrors are now common tools in quantum optics labo- 
ratories HIH. The Bose-Einstein condensates find their 
natural place in this research field as a source of intense 
coherent matter waves. Recently the (coherent) bouncing 
of a BEC on an atomic mirror has been carefully studied 
H . Pursuing the analogy with pure optical phenomena, a 
four wave mixing experiment using BECs has been suc- 
cessfully realized 0. However, contrary to a standard 
four wave mixing of electromagnetic waves, the creation 
of a fourth wave does not require any interaction with 
an external nonlinear medium, but results directly from 
the intrinsic inter atomic interactions. This fact marks a 
dramatic difference between the behavior of matter waves 
and light. 

For future applications of those atomic coherent mat- 
ter waves, it is of great importance to study in detail 
the role of the interactions in the dynamics of the con- 
densate. Such studies are of crucial interest especially 
when using atomic waveguides and atomic resonators. 
The main element of such atom optics tools is a far off 
resonant hollow laser beam which can be loaded with a 



Bose-Einstein condensate from an atomic trap ||. For 
a sufficiently large detuning, spontaneous emission may 
be neglected leading to a coherent manipulation of the 
motion of the atoms. If used alone, this doughnut laser 
beam may be viewed as a matter wave guide, see Fig. 1. 
If two other blue detuned lasers are added transversally, 
an atomic de Broglie cavity can be formed. 



FIG. 
trap. 




1. A schematic view of the condensate in the optical 



In a previous paper || we studied the dynamics of a 
condensate inside a simple resonator (modelled by a one 
dimensional square box) and focused particularly on the 
influence of the inter atomic interactions on the behavior 
of the system. We showed that starting from a localized 
wave packet and by increasing the nonlinearities we pass 
from a quasi linear regime, where phenomena like non- 
linear revivals and fractional revivals are expected Jl0| , 
to a regime where stochastization of the movement takes 
place and where the nonlinear revivals of the wave packet 
are replaced by revivals of regular dynamics. 

In this paper we are interested in finding the same kind 
of phenomena but using a more realistic trapping poten- 
tial. In particular the rigid boundary conditions corre- 
sponding to the square box potential have to be modified, 
since this kind of discontinuities cannot be experimen- 
tally realized. Moreover, in order to probe the dynamics 
of the system, some experimentally accessible tools, more 
adequate than those used in ||, have to be introduced. 
It must be stressed that the (nonlinear) revivals investi- 
gated in this paper concern a mesoscopic object following 
a nonlinear evolution equation. So far the revivals of a 
wave packet have been studied 1 1 , 1^1, and observed only 
for microscopic and linear systems p3 17 . On the other 
hand, the irreversible behavior which may appear in the 
dynamics, concerns the evolution of an isolated system 
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and is only due to the nonlinear character of the evolu- 
tion; it is not caused by any damping originating from 
some interactions with an external reservoir. 

The plan of the paper is the following: in section II, 
after describing briefly the model which we use, we intro- 
duce some experimentally accessible observables, namely 
the variance and the entropy associated with the posi- 
tion, which will be used to check the dynamical behavior 
of the system. We will show how to recover the results 
of Ref. H (obtained for a box potential) in the cases 
most interesting for experiments. In section III we use 
the observables of section II in the case of a more re- 
alistic potential. We discuss how the different regimes 
of the condensate behavior is reflected in the behavior 
of the variance and the entropy. In section IV, in order 
to support the validity of the one dimensional model we 
present some results obtained in a 2D simulation. Fi- 
nally, we conclude in Sec. V. 



II. MODEL AND NEW VARIABLES 

We recall here briefly the mathematical model we use 
in order to describe the evolution of a Bose-Einstein con- 
densate in an atomic cavity. We consider the zero tem- 
perature case and use the standard Hartree-Bogoliubov 
approach Q. All the atoms are described by the same 
wave function ^>(r, t), the evolution of which is given by 
the time dependent Gross-Pitaevskii equation 
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7— V 2 + V tra p(r,t) 

2m 

Nu \V(f,t)\ 2 *(f,t). 



(1) 



We denote by uq — Ai:h 2 a sc /m, with a sc being the s-wave 
scattering length, and by N the total number of atoms. 
The core of the atomic resonator we have in mind is a 
Laguerre-Gauss laser mode. The radial confinement due 
to this optical potential is very strong compared to the 
longitudinal one. Typically the internal radius of the 
"cylinder" is of the order of a few micrometers and the 
length L of the "box" may vary from 75 to 100 microm- 
eters. This leads to a radial confinement of the order of 
several hundreds Hz, while the longitudinal scale of fre- 
quency will be of the order of few Hz. Unless we use a 
very large number of atoms, the dynamics in the radial 
direction can be considered as frozen; the system behaves 
effectively as a ID one, relevant dynamics occurs then in 
the longitudinal direction (we will come back to this as- 
sumption in Sec. IV). Therefore we use a wave function 
in a product form <3>(r, t) — ip(z, t)xo{x, y, t), with xo de- 
noting the ground state in the transverse direction. We 
simplify now the trapping potential by using a square 
box potential. The evolution equation for ip(z, t) is then 
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n 2 d 2 ^{z,t) Nup 2 

g z 2 +— IV>(M)I ^(z,t), (2) 



with ip fulfilling rigid boundary conditions, and S de- 
noting the (effective) transversal section that we assume 
to be constant over the total length of the box. This is 
of course an additional approximation, since the dough- 
nut laser beam is always focused in the region where the 
loading of the condensate takes place. Consequently, rig- 
orously speaking we have S — S(z) with S(z) exhibiting 
a minimum in the center, and a maximum at the borders. 
Since the relative variance of S(z) is rather small |19| , we 
may neglect it in the first approximation. 

Depending on the value of the nonlinear parameter, 
the dynamics described by Eq.(2) can attain very differ- 
ent character The relevant nonlinear parameter is 
given by the ratio where Qi n t = Nuo/h is the char- 
acteristic frequency related to the interaction term and 
u>\ is the frequency of the fundamental mode u>\ = ^mh' 1 • 
We investigate the dynamical behavior for the values 
= 3.08 and = 61 which correspond to a typical 
number of atoms from several hundreds to few thousands 
in the first case, and from few thousands to few tens of 
thousands in the second one. Independently of the chosen 
initial state, we observe that for ^j 2 ^ = 3.08 the dynam- 
ics remains regular and quasi-periodic (nonlinear revivals 
of the wave function occur), whereas for = 61 the 
motion is stochastic. In the latter case the wave function 
does not exhibit any strong spatial relocalization, but 
the previously mentioned nonlinear revivals are replaced 
now by revivals of regular dynamics characterized by a 
return of most of the energy in the few initially popu- 
lated modes. These results were obtained in Ref. || by 
checking the appearance of energy transfer between the 
eigenmodes of the system and by looking at the temporal 
spectrum of the density. 

Experimentally the access to the quantities discussed 
in Ref. || is difficult, if not hardly possible. The sim- 
plest quantity to measure experimentally is the local 
density N\tp(z,t)\ 2 . In order to check the appearance 
of revivals of the wave function or revivals of regular 
dynamics the most natural time dependent quantity to 
study is the variance a z of the mean position of the con- 
densate cloud, defined as a 2 = ((z — {z)^) 2 )^ where 
= / dz...\il>(z, t)\ 2 . If this function exhibits a re- 
turn to its initial value we will be able to conclude about 
the appearance of nonlinear revivals of the wave func- 
tion, or equivalently to some spatial relocalization to its 
initial position. Although this quantity is adequate in 
order to probe the revivals, it is insufficient to conclude 
without ambiguity that a complete derealization of the 
wave function has occurred. Such derealization may cor- 
respond either to a relaxation into another regular state 
with a density peaked only at the cavity ends, or to ir- 
reversible dynamics leading to a more or less uniform 
density distribution. We need therefore another dynam- 
ical quantity which would be more sensitive to various 
kinds of spatial derealization and relocalization of the 
wave packet. The entropy S z associated with the posi- 
tion of the system fulfills such a requirement EG]. This 
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statistical entropy is defined as 



S,(t) = -52P n (t)l0gP n (t), 



with P n being the probability of presence in the interval 

Zn , Zn 



'-m ~r dz] contained in the interval [0, L], 

Pn{t) = \lP(z n ,t)\ 2 dz. 



(3) 



This function reaches its maximum value SJ lax for a uni- 
form probability distribution corresponding to the case 
of a spatially uniform system. On the other hand any 
spatial localization of the wave function will induce a 
strong decrease of S z compared to the maximum value. 
As usual it is better not to use S z but a normalized en- 
tropy r\ z defined as 
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Experimentally, both a z and S z correspond to mea- 
surable quantities. The only restriction is coming from 
the (inevitably) limited spatial resolution in a given ex- 
periment. But this limitation affects only the dynamics 
associated with high frequency modes and the interest- 
ing phenomena (revivals of the wave packet or revivals of 
regular dynamics) are in fact associated to the low modes 
of the system. A typical spatial resolution can be of the 
order 3/im which, in a box of 100/xm limits the resolution 
of the modes higher than the 30th. 

We simulate the evolution of a wave packet initially 
centered in the box with no initial velocity and with an 
initial width of L/5. This corresponds to the first case 
investigated in [^) where mainly the modes 1 and 3 of 
the box potential are populated. We use the usual split 
operator method to solve Eq.(2) ^J. We check the con- 
servation of the norm of the wave function and the con- 
servation of the total energy. The former was conserved 
better than 10 _4 %, and the latter has shown fluctuations 
smaller than 1%. In the figures 2a and 2b we present re- 
spectively the time evolution of the variance compared to 
its value for a uniform system a un i and the time evolution 
of r\ z . These plots correspond to a nonlinear coefficient 

^ La± = 3.08. The time is given in units of the fundamen- 
tal ° 

tal period of the box potential T\ = 2tt/lui. We restrict 
ourself to a time interval of one period roughly since for 
the experimental system it would correspond to a time of 
the order of 18s which is beyond, or in the best case close 
to the usual lifetime of the condensate in this situation. 
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FIG. 2. Evolution of (a) a z /a un i and (b) r\ z in the case 
^ai- =3.08 and for an initially centered wave packet. 



The evolution of a z shows clearly some quasi-periodic 
revivals at times close to those expected in the linear 
case with this initial state. The nonlinear character of 
the evolution appears through the position of this re- 
vival. The behavior of r\ z confirms the regular spatial 
relocalization of the wave function. This function ex- 
hibits clear returns to its maximal value 1. Notice that 
the fractional revivals corresponding to a localization of 
the atoms near the boundaries are not at all exhibited 
by a z , but appear very clearly in the time dependence of 
rj z , as expected. They are represented by the small peaks 
preceding a quasi-revival of the system. The occurrence 
of quasi periodic nonlinear revivals of a mesoscopic wave 
packet is thus confirmed. From the experimental point 
of view the bad point is of course the time scale of these 
phenomena. The first nonlinear revival is shown at a time 
close to O.lTi (~ 1.8s), which is already rather long. We 
will see in section III that in a more realistic potential 
than a box like potential the revival time is significantly 
reduced. 

In the figures 3a and 3b we present the time evolu- 
tion of the same quantities but for a nonlinear parameter 
= 61. 
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FIG. 3. Evolution of (a) a z /a un i and (b) rj z in the case 
Q±at. =61 and for an initially centered wave packet. 



The variance does not exhibit any clear return to its 
initial value. On the contrary it oscillates around the 
value a-uni ■ The absence of strong spatial relocalizations 
of the wave function is confirmed by the behavior of r\ z 
which stays far from 1 and oscillates around a low value. 
Notice that this value is not zero, which means that the 
system has effectively relaxed but is still following some 
dynamics. The change from the nonlinear revivals of tp 
to the revivals of regular dynamics is however visible. 
The occurrence of regular dynamics during some time 
intervals is indeed clearly exhibited by the variance. At 
these times where the energy of the system comes back 
to the low modes of the box potential the variations of 
a z have a significant lower rate as shown in Fig. 3a. This 
kind of behavior appears also in the variations of r\ z but 
is not as spectacular as for a z . 

At this point the variance and the entropy associated 
with the position appear to be good physical quantities 
to probe the dynamical behavior of the system and for 
showing its particularities (revivals of the wave packet 
or of regular dynamics). However, the time scale needed 
for observation of these phenomena is too long from the 
experimental point of view. This fact is mainly due to 
the form of the chosen potential. We investigate in the 
next section how these results are modified by introduc- 
ing more realistic boundary conditions. 



III. STUDIES OF THE DYNAMICS IN A 
REALISTIC TRAPPING POTENTIAL 

Experimentally the optical box will be closed transver- 
sally by using two blue detuned lasers. The trapping 
potential at the boundaries of the box is no more discon- 
tinuous and may be modelled in a good approximation as 
the sum of two identical Gaussians centered at z = and 
z = L. The characteristic parameters of the potential 
are then its height Vi aser and the width of the Gaus- 
sian aiaser ■ The requirement on VJ aser is that it provides 
a sufficiently high barrier in order to prevent the atoms 
from leaving the trap. We take it as a given parameter of 
the problem. The value of ai aser directly depends on the 
transversal focusing of the lasers at z = and z — L. The 
typical value cannot be smaller than a few \xm. Clearly a 
too small aiaser will take us back to the case of the per- 
fect box. It is beneficial to use a large ai aser in order to 
reduce the time scale on which the nonlinear revivals of 
the wave packet or of the regular dynamics appear. On 
the other hand, too large er; aser does not interest us be- 
cause it will correspond to a situation, where the atoms 
are still very confined and the cloud will just oscillate as 
in a slightly modified harmonic trap. In other words, the 
overlapping of the boundary barriers with the initial state 
must remain as small as possible in order to remain in a 
situation similar to a box potential from the dynamical 
point of view. 

The study of the dynamics of the system appears to 
be simpler if one uses the quantities rj z and a z , since it is 
difficult to have access to a large range of different modes 
of the potential (which would be necessary in order to 
compute for instance the statistical entropy associated 
with the energy distribution among the eigenmodes as in 
Ref. ||). We present here the results of the simulations 
for the same case as before, namely an initially centered 
wave packet with a zero mean velocity for the two same 
values of the nonlinear parameter. The width of the wave 
packet is still L/5 and we choose a ai aser = L/12 which 
corresponds in a typical experiment to a focused waist of 
roughly 16/Ltm and the height of the barrier is chosen to 
be of the order a few mK. We show the time variations of 
the same quantities than before in order to compare the 
results with the case of the perfect box. The variance 
&z/Puni is still expressed in units of the uniform case 
which is calculated using the box potential. We keep also 
the same unit of time which will allow us to compare the 
time scale obtained now for the different phenomena we 
want to investigate. In Figs. 4a and 4b we show the 
variations of the variance and the normalized entropy r\ z 
for Ote* = 3.08. 
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FIG. 4. Temporal variation of (a) a z /a un i and (b) i) z 
for an initial centered state without mean velocity where the 
boundary potential has a Gaussian shape and ^j^ 1 = 3.08 

The regular and quasi-periodic behavior of the two 
quantities is striking. The wave packet clearly exhibits 
nonlinear revivals at almost regular times. Although the 
maximal value of a z is as expected smaller than for the 
perfect box, the most important point is that the time 
scale for the nonlinear revivals has been divided by a fac- 
tor 2. This makes the first quasi-revival of the condensate 
to appear before Is. Although it remains long, it is now 
sufficiently shorter than the lifetime of the condensate. 
This time scale can be decreased even more by playing 
with the value of <ri aser as well as with the initial state 
of the system. In particular a spatially more squeezed 
initial wave packet will populate higher excited states at 
the beginning of the evolution leading to a smaller lin- 
ear revival time, and, consequently, to a smaller nonlin- 
ear revival time (with these values of the nonlinearities 
the two regimes are not too much different and the time 
scales are comparable). It is important to note that, as 
before, some fractional revivals corresponding to an ac- 
cumulation of atoms at the boundaries are still visible in 
the variations of r\ z . 

We turn now to the case of strong nonlinearities. The 
simulations of the evolution a z and r\ z for a nonlinear 
parameter = 61 are presented in Figs. 5a and 5b 
keeping the same conventions as before. 
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FIG. 5. Temporal variation of (a) a z /a uni and (b) rj z for an 
initial centered state without mean velocity where the bound- 
ary potential has a Gaussian shape and ^ =61 

The behavior is similar to those of Fig. 3. The variance 
does not come back to its initial value and r\ z stays far 
from 1 which proves that the wave function does not 
exhibit any spatial relocalization. The variations of a z 
still exhibit the revivals of regular dynamics. Although 
the time scale on which they appear has been reduced it 
is still of the order of O.lTi which makes them difficult 
to observe experimentally. 

IV. ANALYSIS OF THE TWO DIMENSIONAL 
PROBLEM. 

A. Analysis of the model. 

In the previous sections we have studies the simplified 
ID dynamics arguing that the radial confinement due to 
the doughnut laser mode was much stronger than the 
longitudinal one. In consequence the structure of the 
condensate in the radial direction will not be modified 
by the nonlinearities. Although this hypothesis is rea- 
sonable, it has to be checked at the level of the physical 
quantities we are using to describe the behavior of the 
system. The ratio of the typical energies associated with 
the radial degrees of freedom and the longitudinal ones is 
of the order of L 2 /R 2 ~ 100 where R denotes the internal 
mean radius of the Laguerre-Gauss laser mode and L the 
length in the longitudinal direction. Qualitatively, for 
small nonlinearities, when the dynamics is already quasi- 
periodic in the longitudinal direction, there are no serious 
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worries concerning the influence of the other degrees of 
freedom. On the other hand, for stronger nonlinearities 
corresponding to of the same order of magnitude as 
L 2 /R 2 , a strong coupling between the different degrees of 
freedom associated with different directions may appear, 
and it is not sure that we can get the same behavior as 
before for the variance and the entropy associated with 
the longitudinal direction. However, we must keep in 
mind that this simple energetic estimate does not take 
into account the population of the modes which are im- 
portant parameters in nonlinear problems. The deriva- 
tion of more accurate criteria is then necessary. We will 
propose some of them in the next paragraph. 

We turn now to the discussion of the inclusion of an 
additional dimension. For this purpose we consider the 
evolution of a condensate in a two-dimensional box. We 
denote the directions by z and y and by L z and L y the 
associated lengths. We take L z = 4L y where the energy 
scale in the y direction is 16 times bigger than the one of 
the z direction and will mimic our radial confinement in 
the real system. Of course as we are working in Cartesian 
coordinates the problem is simplified compared to the 
cylindrical one but it should give us a good representation 
of the real system. Moreover the typical ratio of the 
energies is one order of magnitude smaller than the one 
of our system but it will still provide the same kind of 
behavior and allows us to test values of nonlinearities 
sufficiently small to have a good spatial resolution in our 
numerical grid. 

The starting point is still the time dependent Gross- 
Pitaevskii equation (Q). Assuming a very strong con- 
finement in the x direction, we can neglect the degrees 
of freedom associated with the x direction by assuming 
a wave function of the form \&(r, t) = ip(z,y,t)£ (x,t) 
where £o is the fundamental state of the x direction. The 
evolution equation of ip is then given by 
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— V 2 ^(z,y,t) 
2m 

\^(z,y,t)\ 2 ip(z,y,t), 



(5) 



where L x denotes the length in the x direction. 

As in the one dimensional case in 1J we can derive 
some criteria to estimate the kind of dynamics the sys- 
tem will have. We repeat here briefly these ideas. We 
decompose ip using the basis of the eigenmodes of the two 
dimensional box, ip(z,y,t) = Y,n, m c nm.{t)4'n{z)xm{y) 
with cj> n {z) = yX sin (-L£) and Xm ( y ) = y|T s i n (I!^/). 

Inserting this form into the equation (^) we obtain a 
set of coupled differential equations for the coefficients 
c nm . They can be derived as the Hamiltonian equa- 
tions of the motion for the complex degree of freedom 
{c nm ,ihc^ m }. It is however more interesting to work 
in action-angle variables {/ n m,^nm} defined such that 

c nm = \ ^MQ-e ie " m . The Hamiltonian leading to the dif- 



ferential equations of motion for the c nm can the be writ- 
ten in action-angle representation as 



n m 

+ 1 



^ ^ ^nm(^nlml^n2m2-^n3m3-^n4m4) 



/2 



~i {@nlml-\-6n2rn2—8n3rn3—8ri4m<i) 



(6) 



where q = -t-ttt and where 



^n.rri. 



L Z f-Ly 
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Here ft and m correspond to (711,77,2,77,3,77,4) and 
(mi, m2, 777,3, 7714). We have also introduced the frequen- 
cies uj nm corresponding to the mode labelled (71,777). 
We have of course huj nm = " 2 ^^ 2 + "f^lf = {n 2 + 

I6m 2 )hu>i z where we introduced the fundamental fre- 
quency LOi z in the z direction. 

Compared to the equations derived in j^] we have just 
a doubling of the indices due to the inclusion of a second 
dimension. The derivation of the criteria which indicate 
the borderline between two different dynamical behaviors 
is then straightforward. The first criterium indicates for 
what conditions one may expect that the dynamics of the 
systems will remain close to the linear one. In order to 
derive this criterium we have to compare the frequencies 
<-<Jnm of the linear problem with their first order correction 
given by the nonlinearities |]] . In this case we have 
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Ti(n 2 + 16m 2 ) 



(7) 



where the natural frequency of the nonlinearities is de- 
noted by flint = ^p- As the I„ m 's are normalized to 
h the high modes of the box will satisfy this inequality 
more easily than the low ones for a given value of the 
nonlinear parameter Qi^. 

The Chirikov criterium gg^] on the other hand, gives 
us an indication about the conditions for which the dy- 
namics may stay regular. This condition follows from 
the resonance overlap criterion which expresses the con- 
dition that separatrices associated with two consecutive 
resonances are touching each other. We compare the dif- 
ference between two consecutive bare frequencies with 
the frequency associated with the first order correction. 
In our case it leads to 



, 16 uu w 
< — — — n(n 

9 ilint. 



16m 



2 1 



(8) 



As the previous inequality, this one is also more easily 
fulfilled for the highly excited modes, which consequently 
may preserve a regular motion easier than the low ones 
for an increasing value of the nonlinearities. Clearly for a 
given nonlinear parameter which is at best of the order of 
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the typical ratio of the two energy scale (say less than 30), 
only the modes with very low m may not satisfy these 
inequalities. As the initial state we choose is a centered 
wave packet of widths L z /10 in the z direction and L y /10 
in the y direction without any mean velocity, we initially 
populate the modes (n, m) with n,m =1 or 3 essentially 
This choice is coherent with what has been done in the 
one dimensional problem and is one of the most favorable 
in order to see the appearance of characteristic nonlinear 
dynamics. 

We consider three cases for which the ratio takes 
the values: (a) 4, (b) 10, (c) 101. Experimentally this 
would correspond to condensates whose number of atoms 
are: (a) few hundreds of atoms, (b) from few hundreds 
to 2000, (c) from few thousands to 15000. As for the one 
dimensional case, we want to compute the variance and 
the statistical entropy associated with a given direction. 
In the following paragraph we give the definition of a z 
and S z . Exchanging z «-> y in the following expression 
gives the expressions of a y and S y . In the case of an 
initially centered wave packet, the preservation of the 
symmetry of the wave function simplifies a little bit the 
expressions. We obtain 



1.25 



!(*) 



dz z 2 Tl(z,t) 



(9) 



where 11(2, t) is the density of probability of presence at 
z at the time t, H(z,t) = f Q y dy \tp(z,y,t)\ 2 . For the 
entropy we use the probability distribution P(z n ,t) = 
U(z n ,t)dz, giving S z (t) = - £) fc P(z k , t)logP(z k , t), 
where the dz is just the numerical step used to define 
the grid. As before it is more convenient to work with a 
normalized entropy r\ z whose definition is the same as (Q) 
and with a variance compared to its value in the uniform 
case. 



B. Numerical results 

We have numerically solved the two dimensional time 
dependent Gross-Pitaevskii equation (|B|) using an explicit 
method |2jJ . We checked the conservation of the norm of 
the wave function (better than 10 -4 ) and of the energy 
(better than 1%). In the Figs. 6 and 7 we present the 
time variations of the variance and the normalized en- 
tropy associated with the different directions for the case 
(a). The time is given in unit of T\ = 2tt/loi z . This is 
the reason why the dynamics in the y direction is pre- 
sented on a time scale ten times shorter than for the z 
direction. As in the one dimensional case, it is however 
still long from an experimental point of view. 
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FIG. 6. (a) Temporal variations of a z /a un i Z ; (b) Temporal 
variations of a y /u uniy for an initial centered state without 
mean velocity for that. — 4 
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FIG. 7. time variations of the normalized entropy associ- 
ated with the position for an initially centered wave packet, 
(a) v „ (b) Vy for = 4 
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The dynamics is regular as we may expect it for this 
small value of nonlinearitics. For times shorter than 
0.057i, the behavior associated with the y direction is 
however closer to a linear one than in the z direction. 
This appears through the returns of the normalized en- 
tropy to the value 1 as it will be the case for a linear 
behavior. This means that the m mode distribution did 
not change during the evolution on the time scale pre- 
sented here. For longer times a decrease in the precision 
of the revivals appears, but it remains rather small. In- 
deed, it is only visible on the scale of the normalized 
entropy although the variance did not exhibit a notice- 
able change. In the z direction quasi-periodic revivals of 
the wave packet in this direction are shown. They corre- 
spond exactly to the ones seen in Fig. 2. We can then 
conclude that for the smallest value of the nonlinearitics 
considered in the one dimensional case, the effects asso- 
ciated with the radial degree of freedom will not prevent 
the appearance of nonlinear revivals of the wave function, 
and in fact do not play a relevant role. 

We turn now to the case of intermediate nonlineari- 
ties corresponding to the case (b). In Figs. 8 and 9 we 
present the variations of the variances and the normalized 
entropies. 
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FIG. 8. (a) Temporal variations of a z /a U niz; (b) Temporal 
variations of a y /a uniy for an initial centered state without 
mean velocity for Q±m, = \q 
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FIG. 9. time variations of the normalized entropy associ- 
ated with the position for an initially centered wave packet, 
(a) Vz , (b) Vy for ^nt = 10 



The dynamics associated with the two directions are 
completely different. In the z direction the quasi-periodic 
spatial relocalizations of the wave function have vanished. 
As in the one dimensional case, they have been replaced 
by revivals of regular dynamics as can be seen in the in- 
set for the variance a z . A stochastic dynamics has taken 
place for this direction. On the contrary the dynamics in 
the y direction has remained regular and quasi-periodic, 
although compared to the previous case the nonlineari- 
ties begin to influence it stronger. Consequently when 
the nonlinear parameter is comparable to the ratio 
between the energies in z and y direction, revivals of reg- 
ular dynamics in the direction of the small energy scale 
are expected. The predictions done in the approximated 
one dimensional case are then valid for the typical values 
of the nonlinearities we have considered (we point out 
that in this case we have investigated the behavior for 
— 61 where the ratio of the scale of energies is of 
the order of 100). 

If we increase more strongly the nonlinearities the cou- 
pling between the modes can no more be neglected for the 
y direction. This is shown in the Figs. 10 and 11, where 
we present the variance and the normalized entropy for 
the case (c). 
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FIG. 10. (a) Temporal variations ol cr z /a un i Z ; (b) Tempo- 
ral variations of o- y /a un i y for an initial centered state without 
mean velocity for Qiat- — \q\ 
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FIG. 11. Time variations of the normalized entropy associ- 
ated with the position for an initially centered wave packet, 
(a) Vz , (b) Vy for %*f = 101 



The very fast decrease of the two normalized entropies 
indicates a strong relaxation of the system towards a 
nearly uniform state. The variance in the z direction 
reaches almost a steady state. It oscillates around its 
value corresponding to the uniform case and the revivals 
of regular dynamics of the previous case are still visible, 
but vanish very quickly. In the y direction the situation 
is similar. The previous revivals of the wave packet in 
this direction have been replaced by revivals of regular 
dynamics, whose amplitudes decrease quickly with time. 
The dynamics have developed stochasticity as may be 
expected from the application of the Chirikov criterium 
(Q) for this value of the nonlinear parameter. 



V. CONCLUSIONS 

In this paper we have studied the nonlinear dynamics 
of a condensate trapped in an optical box like potential. 
The scheme we have in mind is a condensate located 
in the center of a box-like potential where we monitor 
the time evolution of the wave packet by looking at the 
associated entropy and variance in order to detect the 
main properties of the dynamics. The behavior of these 
functions allows to confirm the predictions done previ- 
ously in ||, namely that for small nonlinearities (corre- 
sponding to situations with few thousands of atoms at 
most) the system will exhibit nonlinear revivals through 
quasi-periodic spatial relocalizations. On the other hand, 
for stronger nonlinearities (corresponding typically to a 
number of atoms of few tens of thousands ) the system 
will follow an effective relaxation dynamics, for which the 
wave function revivals are replaced by revivals of regular 
dynamics. The observed revival times in a perfect box 
potential are shown to be rather too long to be exper- 
imentally measurable. Fortunately using a more realis- 
tic trapping potential with a pair of blue detuned lasers 
as "end caps" for the doughnut mode laser configura- 
tion, there is a substantial reduction of the revival times, 
which can then be shorter than the typical condensate 
life time. It is important to note that the present study 
demonstrates the possibility of observation of the revivals 
of a mesoscopic object, which is interesting from an ex- 
perimental point of view. It shows, however, also irre- 
versibility in a closed system which is interesting from a 
more fundamental point of view. 

In the numerical treatment we have first assumed all 
the dynamics was taking place in one dimension. In other 
words the transversal dynamics of the condensate was 
frozen. This assumption, based on an effective strong 
radial confinement, is valid as long as the (nonlinear) 
coupling is not too strong. In order to justify the one di- 
mensional treatment and to study the break-down of it, 
we have performed a two dimensional simulation using 
a perfect box potential. By choosing a trap aspect ratio 
significantly lesser than one, we have demonstrated that 
the coupling between the two dimensions is present only 







at large values of the parameter f2,- n j /w\ z describing the 
nonlinear interaction. Hence, the one dimensional study 
provides indeed a good approach for the physical parame- 
ters used here. The coupling between the two dimensions 
is mainly governed by the different energy scales of the 
eigenmodes in z and y direction due to the aspect ratio. 

Finally, we would like to point out that very recent ex- 
perimental progress of the Hannover group allows to ex- 
pect that the observation of the condensate wave packet 
and transition to irregularity will be feasible within the 
times of the order of Is [ pj[ . 
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